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I. INTRODUCTION 



In order to get a complete information about quantum mechanical systems, one should 
study bound and scattering states in the presence of an external potential. Therefore, 
scattering problem has become interesting topic in relativistic/non-relativistic quantum me- 
chanics. The problem has been studied within the framework of group theoretical approach 
[1] and investigated for the well-known potentials by applying different methods [2-18]. In 
the case of non-relativistic scattering problem, it has been showed that transmission and 
reflection coefficients take 1 and 0, respectively, as external potential has well-behaved at 
infinity for the zero energy limit [10-12]. However, reflection coefficient goes to zero while 
transmission coefficient goes to unity in the zero energy limit when external potential sup- 
ports a half-bound state. This situation was called as transmission resonance by Bohm [13]. 
This phenomenon has been recently extended to the Dirac particle in Rcfs. [3, 14, 15]. In 
the zero-momentum limit, Dombey et. al displayed that the bound-state energy eigenvalue 
obtained for the Dirac particle in the presence of the Woods-Saxon potential is related to the 
transmission resonance appearing for a Dirac particle scattered by a potential well. Recently, 
Villalba et. al [6] have showed that relation between the bound-state energy eigenvalues and 
transmission resonances in view of the Woods-Saxon potential for the Klein-Gordon particle 
are the same as obtained for the Dirac particle [15]. 

On the other hand, solutions of the wave equations have, recently, become interesting 
in the view of position-dependent mass (PDM) formalism. Extensive applications of this 
formalism have been done in different areas of physics such as condensed matter physics 
and material science such as electronic properties of semiconductors [19], quantum dots 
[20], quantum liquids etc. [21-25]. In recent years, the scattering problem has been ex- 
tended to the case where the mass depends on spatially coordinate [26-28]. Alhaidari has, 
recently, investigated solution of the Dirac equation in view of position-dependent mass for 
the Coulomb field [24]. In Ref. [27], the authors have studied the relativistic scattering in 
the Dirac equation by using the J-matrix method for the position-dependent mass. Panella 
et. al [28] have obtained a new exact solution of the effective-mass Dirac equation for the 
Woods-Saxon potential. The approximate solution of the Dirac equation with PDM for the 
generalized Hulthen potential has been obtained by Peng et. al [29]. Jia and co-workers have 
extended PT-symmetric quantum mechanics for the Dirac theory to the PDM formalism 
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[30, 31]. 

In this work, we intend to solve the effective-mass Dirac equation for the Woods-Saxon 
potential and investigate the scattering and bound state solutions. We also study the trans- 
mission resonances and give some results in the case of the low momentum limit. 

The organization of this work is as follows. In Section 11, we give the one-dimensional 
Dirac equation for the case of position-dependent mass and obtain a Schrodinger-like equa- 
tion. In Section III, we calculate the transmission and reflection coefficients by analyzing the 
behavior of the wave functions at x — > ^oo. To compare our results, we also calculate the 
same coefficients for the case of constant mass and give our results. In Section IV, we study 
the bound state problem for the effective-mass Dirac equation and give the transmission 
resonance and bound state equations for the low momentum limit. Conclusions are given in 
Section V. 



II. DIRAC EQUATION WITH POSITION DEPENDENT MASS 

The relativistic free-particle Dirac equation (h — c—1) is written as [3, 15] 

[iji'd^ - m{x)] iP{x) = , (1) 

where we assume that the mass of the Dirac particle depends only on one spatially coordinate 
X. Under the effect of an external potential V{x) and taking the gamma matrices ^y^ and 70 
as the Pauli matrices ia^ and a^, respectively, the Dirac equation in one-dimension becomes 

which gives the two following couple differential equations 

d(pi{x) 




dx 
d(fi2{x) 



-[E-V{x) + m{x)]ip2{x), (3) 
[E - V{x) - m{x)] (pi{x) . (4) 



dx 

The solutions can be more easily obtained by using a two-component approach introduced 
by Fliigge as [32] 

(f){x) = (pi{x) + i(p2{x), (5) 
X(x) = (pi{x) -i(p2{x), (6) 



which leads 

= i[E - V{x)](p{x) - im{x)x{x), (7) 

= -i[E - V{x)]x{x) + im{x)(j){x). (8) 

Ehminating in Eq. (7) and inserting into Eq. (8) and following the similar procedure 
for (f){x), we obtain two uncoupled second-order differential equations for (f){x) and x(x), 
respectively 

d?(j){x) dm{x) / dx d(j){x) , 2/ \ .dV{x) 



+ S^[E-V{x)]'' -m^{x)+^- 



dx'^ m(x) dx I dx 

n{x)/ 
m{x) 



+ z[E - L^(^) ^ , (9) 



d^x{x) dm{x)/dxdx{x) ^ f _ y (^^^^^2 _ ^2 (^^^^ _ -dV (x) 



and 

a 

dx'^ m{x) dx ' \^ ^ ^ ' dx 

-i[E-V{x)\'^^^, 

m{x) ) 

III. MASS FUNCTION AND SCATTERING STATE SOLUTIONS 

We assume that the mass of the Dirac particle depends on spatially coordinate giving as 

m{x) ^ mo + mif{x) , (11) 

where the function of x as f{x) — 1/(1 + e"^'^'"^-*). The parameter mo will correspond to 
the rest mass of the particle and mi is a real, positive, small parameter. The mass form 
provides us to obtain the analytical results for the reflection and transmission coefficients 
and also the bound state solutions for the case of position-dependent mass and to analyze 
the results for the case of constant mass. On the other hand, from Eq. (11), it is easy to 
see that the ratio of the derivative of the mass to the mass is proportional with the mass 
parameter mi. So we ignore the terms that contain the derivative of the mass in Eqs. (9) 
and (10) for the case of mi — > [29]. Under this assumption, Eqs. (9) and (10) become 

^''^(^) + ![E- V{x)]' - m\x) + z^^\<p{x) = , (12) 



dx'^ I dx 
dx^ 



+ i^[E-V{x)f -m\x)-i^^^^x{x)^0. (13) 



We search the scattering states of the Dirac equation for the WS potential barrier [3] 

V{x) = Vofix) , (14) 

where the function f{x) is defined in Eq. (11) with Vo,a and L are real parameters. This 
potential is one of the most important potential models in quantum mechanics and has 
a main role, as an internuclear potential, in the coupled-channels calculations within the 
heavy- ion physics [33]. The nuclear optical- model potential including the Woods-Saxon 
potential is used to analyze the elastic scattering problem of nucleons and heavy particles 
[34]. It is worth to say that we deal with a potential form for aL ^ 1. In this case, the 
potential form now closely becomes a rectangular barrier. 



A. Solutions for x <0 

Using Eqs. (11) and (14) and defining a new variable y = {1 + e-"(^+^))-S Eq. (12) 
turns into 

^(1 - + (1 - + ^ {C. - C.. + C^y^} Mv) ^ , (15) 

where 

-ml _ 2VoE + 2momi - iVoa _V^ -m\- iVoa 

— n ; (^2 — 5 'i^3 — o ■ V-I-OJ 

q;^ 

In order to get a hyp ergeometric- type differential equation, we offer a trial function (piiy) = 
y^{l — yY f{y)- Substitution it into Eq. (15) leads to 

y{l-y)'^^ + [l + 2^^-2{^, + v + l)y]'^-{^^ + v + a){^l + v-a+l)f{y) = Q, 

(17) 

with 



(18) 



Eq. (17) has a general solution 

f{y) = Li 2-Fi(/i + + cr, // i/ - (7 -Fl, 1 2//; y) 

+ L2y-''^2i^i(-A* + i^ + (T,-At + i/-(7 + l,l-2//;y), (19) 
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which gives 



My) = W(l-y)%Fi(a,6,c;y) + L2y-'^(l-y)%Fi(ai,6i,ci;y), (20) 

where 

a = ii + v + a ; ai = — + + ex , 
h — — (T-\-l] hi — —/I + u — a + 1 , 

c = 1 + 2// ; ci = 1 - 2//. (21) 



B. Solutions for a; > 

In this case, changing the variable z = 1/(1 + e°'^^~^^) and inserting Eqs. (11) and (14) 
into Eq. (12), we get 

41 - z)"^ + (l-2z)'^ + {C. - C,z + C,z^} MA = . (22) 

with 

_E'^ - ml _2E + 2momi + iVoa _V^ -ml + iVoa 

(^4 — 5 — ; (^5 — 5 ; ^6 — n ■ i^-jj 

a"^ 

Defining a wave function of the form (j)R{z) — z'^{l — z)''h{z) in Eq. (22) gives a 
hypergeometric-type equation [35] 

z(l-^)^^^ + [l + 2T-2(r + 7 + l)^]^^-(T + 7 + 5)(T + 7-5 + l)/i(^) = 0, 

(24) 



where we have obtained as t — n ,^ — v and <^ = | + y — |) + ^ ■ Eq. (24) has a 
solution in terms of hypergeometric functions [35] 

h{z) = Ri2Fi{n + v + 5,n + v-5 + l,l + 2ix]z) 

+ R2Z-^^' 2Fl{-^l + v + 5,-il + v -6+1,1- 2/x; z), (25) 

which gives the following complete solution for x > 

(t>R{z) = RiZ>'{l-zy2Fi{a^M.c^\z) + R2Z-''{l-zy2Fi{a2M.C2\z), (26) 
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where 



a2 = -)U + z/ + 5; ^ iJ, + u + 6 , 

62 = —/jl + u — S + I ; bs^ii + u — S + 1, 

C2 = 1 - 2/i ; C3 = 1 + 2/i . (27) 

li X ^ 00, z ^ 0, then (1 - 2;)'' ^ 1 and 2;^* ^ g-a;i(x-L) 'pj^^g^ obtain the following 
right solution in this limit 

(f)R{x) R^e-'^^''-^^ + R2e'^^''-^\ (28) 

where we have used the following property of the hypergeometric functions: 
2-^1(^1)^2,^3;^) 1- III order to get a plane wave coming from the left to the right, 

we set i?i = 0. Consequently, the right solution becomes 

(t>R{z) ^ R-2Z-^{1 - zY 2Fi{a2M,c-2]z). (29) 



C. Reflection and Transmission Coefficients 

Let us now study the behavior of the wave functions (/>l(|/) and (I)r{z) at infinity to obtain 
the reflection and transmission coefficients. In the limit )■— oo,i/— )-0, (1— y)'^— > 1 and 
yii _^ ^a^l(x+L) ^YiQ limit a; — )■ 00, ^ — )■ 0, (1 — z)" 1 and z'^ ^-afiix-L) 

(with 2-^1(^1)^2)^3)^) 1)) we have the wave functions, respectively, 

0l(,t) > Lie''=(^+^) + L2e-*'=(^+^) , (30) 

X — > —00 

M^) > i?2e''("-'') . (31) 

a; — >■ 00 

In order to get the electrical current density for the one- dimensional Dirac equation 
defined by 

J = lm^)f-\x{x)\'], (32) 
we need to insert Eqs. (30) and (31) into Eq. (7) which gives Xl{x) and Xi?(a;), respectively, 

Xl(x) = [{E - A;)Lie^'=(^+^) + (E + A;)L2e-^^(^+^)] , (33) 

Xn{x) = (^) R2e^^^^^^\ (34) 
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The current in Eq. (32) can be written as Jl — jm — jrefi in the hmit x — > — oo where jin is 
the incident and jrefi is the reflected current. Similarly as x — > oo the current is Jr — j trans 
where jtrans is the transmitted current. Inserting Eqs. (30), (31), (33) and (34) into Eq. 
(32), we find the refiection and transmission coefiicients, respectively, as 

r=|f^. (36) 

We can find more explicit expressions for the above coefficients by using the continuity 
conditio of the wave function at x — 0. For the limit x — > 0, we have y — > 1 and 
2-^1(^1)6)6;*) 1) the wave function 0l(x) (oL > 0) 

(I>l{x) > iLiSie-"""^ + 12826-'"'^] e-'*^^ + [Li^gc"'^-^ + 12846'"'^] e'*^^ , (37) 

where we have used the following identity of the hypergeometric functions [35] 

2-^1(6)6)6;*) = (6) — |ll 2Fi(^i,{2)6 + 6 - 6; 1 - *) 

i i?3 - ?lji (,?3 - ?2j 

+ (1 - t)g3-6-6 r(6)r(6 + 6 - 6) ^^^^^^ _ _ - 6 - 6 + 1; i - *) • 

i Uiji i6j 

(38) 



The abbreviations in Eq. (37) are 



_ r(c)r(c-a-6) _ r(ci)r(ci - oi - 61) 



r(c — a)r(c — 6) ' r(ci — ai)r(ci — 61) ' 
_ r(c)r(a + 6 - c) _ r(ci)r(ai + 61 -ci) 
- r(a)r(6) ' fK)r(M ' ^ ^ 

We have 2; — >■ 1 for the same limit x — >■ 0, so we write (f)R{x) as 

M^) > R2 [^5e"^("-^) + S'ee-"^^^-^)] , (40) 

where 

_ r(c2)r(c2 - a2 - ^2) ^ _ r(c2)r(a2 + &2 - ca) , . 

' r(c2 - a2)r(c2 - 62) ' ' r(a2)r(62) ' ^ ^ 



Finally, from matching the wave functions in Eqs. (37) and (40), we obtain 

L2 e "^""^^SiS^ — S^Sq 



(42) 
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and 



(43) 



In Fig. 1, it is seen that the transmission and reflection coefficients oscillate between the 
values zero and one and satisfy the condition i? + T = 1 for the case of constant and PDM as 
well. The oscillations appear in the range of mo < -B < 2mo while T — {R — 1) as E < mo 
and T = 1 (i? = 0) as > 2mo. The effect of the PDM is just to shift the picks to the left. 
We present the variation of the transmission coefficient with respect to the parameter Vq in 
the case of PDM in Fig. 2. We also plot the same variation for the constant mass. The 
transmission coefficient goes to zero with increasingly high of potential barrier and exactly 
zero for the values in the range 0.4 < Vq < 1.2. The coefficient T starts to oscillate and 
does not take zero- value for Vq > 1.2 where the upper value of the oscillation for T is one. 
Fig. 2 shows also the effect of PDM on the dependence of T on potential parameter Vq and 
that this effect is very weak. Fig. 3 shows the dependence of the transmission coefficient on 
the potential parameters a (left plot) and L (right plot) in the case of PDM in view of the 
varying particle energy. In both of the plots, T oscillates between the values zero and one 
as in the case of constant mass. It is seen that the frequency of the oscillations increases 
while the parameter L increases. In Fig. 4, we show the effects of the potential parameters 
a (left plot) and L (right plot) on the variation of the transmission coefficient with varying 
potential parameter Vq. In these figures, coefficient of transmission is exactly zero within 
the range of mo < Vq < 3mo and continue to oscillate out of this range. 

One interesting point is the so-called "transmission resonances" appearing especially in 
relativistic domain [3, 6, 14, 15]. Within our present formalism, considering Eq. (42), the 
transmission resonances {R — 0, T — 1) occur when 



From Figs. 1 and 2, we see that the Dirac particle has transmission resonances in the 
case of PDM like the case of constant mass. In both of figures, we observe that the effect of 
mass depending on coordinate is to shift the picks to the left. Fig. 3 shows that, as in the 
case of constant mass, transmission resonances appear. From the left plot of Fig. 3, one can 
observe that the width of the resonance peaks decreases and the number of the transmission 
resonances remains the same while the parameter a decreases. In addition, it is seen from 



(44) 
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the left plot of the Fig. 3 that the first resonance peak appears at smaller values of the Dirac 
particle's energy in the presence of the small value of the parameter a. Prom the right panel 
of the Fig. 3, one observes that number of the resonance peaks decreases while the width of 
the resonance peaks increases with decreasing the parameter L. From Fig. 4, we see that 
transmission resonances could be observed with varying Vq in the case of PDM and width 
of the resonance peaks decreases while the parameter a becomes smaller but it increases as 
the parameter L decreases. 



IV. BOUND STATE SOLUTIONS 

We tend to find the bound states for the Woods-Saxon potential well which means Vq — > 
-1/oinEq. (14). 

A. Solutions for a; < 

In order to get a complete solution for this region, we use a new variable y = [1 + 
g-a(a:+L)j-i ^^2) and take into account Vq — >■ —Vq, we have 

.(1 - + (1 - 2*)^ + RT^ {'^l - ^'^^ ^ ^3.^} My) = . (45) 

where 

, _ _ ^2 ^ -2VqE + 2momi + iVoa , _ -ml + iVoa 

— 2 ' ^2 — 2 ' ^3 — 2 ■ V*"'' 

a CK 

C{=Ci C!y >C2 Ck >C3 

Vo^-Vo Vb->-Vb 

Taking a trial wave function as (piiu) = y^' — vY' g{y) and inserting it into Eq. (45) we 
obtain 

y(l - y)'^ + [1 + 2/.' - 2{„' + + - (/^' + ^' + '^W + i^'-<7' + l)g{y) = , 

(47) 

with 



(48) 
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The solution of Eq. (47) is written in terms of the hypergeometric type functions [35] 



My) = L,y^\l-yr\F^{a',h\d-y) + L,y-^'{l-yY 2F,{a\X,c!^\y) (49) 



where 



a' — n' + u' + a' ; 
6' = + -cr' + l ; b 
c' = 1 + 2/i' ; 



= —fi + u + a , 
^ -H' + u' -a' + 1, 
= 1 - 2/i' . 



(50) 



B. Solutions for a; > 



Inserting the potential function 



Vo 



1 + ea(x-L) ' 

into Eq. (12) and using the variable z — + e"^^"^)], we get 



dz z{l — z) 



with 



I _ - niQ , _ -2VqE + 2momi - iVoa _Vq -m\- iVoa 



4 ^2 ' 5 

Oi 



C'=Ci 



C' 



(51) 



-(1 - -)^^ + (1 - 2^)^ + TTT^ 1^4 - ^^^^ + C',z'} Mz) - , (52) 



(53) 



Taking a wave function of the form (j)ji{z) — ^'^'(1 — zy'w{z) in Eq. (52) gives a 
hypergeometric-type equation [35] 



z{l - z)^-^ + [1 + 2/x' - 2(/x' + + l)z]^^ - (/.' + + W + u'-5' + l)w{z) - , 



dz^ 



dz 



(54) 



where = | + v ~ |) + ^ ■ The general solution of Eq. (54) is written in terms of 
hypergeometric functions as follow 

w(z) = R3 2Fi{ii' + i/' + 5', At' + i/' - y + 1, 1 + 2/x'; ^) 

+ R^z-''^'' 2Fi{-ii' + + S', -11' + u'-S' + l,l- 2ii'; z), (55) 



11 



and the whole solution for x > is 



M^) = R^^'^'i^ - ^Y' 6'3, 4; z) + R^z->''{1 - zf 2Fi(4, b'^, 4; z) (56) 



where 



4 = 



— /i + Z/+0 ; a3 = /i+z^+o, 
1-2//' ; 4 = 1 + 2//'. 



(57) 



Let us now extract the solutions given in Eqs. (49) and (56) in the limit x — > ±00 to obtain 
the bound state wave function. Because of y — > and 2; — > as well for the limit x — > +00, 
we write 0^ in Eq. (49) and (f)R in Eq. (56), respectively 



My) > W(l - vY 2Fi(a', 6', c'; y) , 



(58) 
(59) 



where we set L4 = i?4 = for obtaining the bound state eigenf unctions. In order to study 
the behavior of the solution at x = 0, we need the property of the hypergeometric functions 
given in Eq. (38). Recalling that for x — )■ 0, z — )■ 1 and (1 — y)"' g-"'^'(a;+^) -while 
(1 — zY' ~ and using Eq. (38), we write the wave functions 



<i>L{x) 



^-av'{x+L) _|_ ^1 ^av' {x+L) 
g/^^au'(x-L) ^ g,^^-au'{x-L) 



(60) 
(61) 



where 



_ T{d)T{d-a'-h') _ _ r(4)r(a'3 + 6^ - 4) 



V{d - a')V{d - h') 



r(a^)r(4) 



r(4)r(4 - 4 - 4) 



^ V{d)T{a' + b'- d) 

T{a')T{h') ' ' r(4-agr(4-4) 



(62) 



Matching the functions given in Eqs. (60) and (61) at x = requiring the continuity of the 
wave function, comparing the coefficients of e^""^'^ and setting the coefficients determinant 
to zero, we obtain the following eigenvalue condition for the Woods-Saxon potential well 



/[a, Fo, mo, mi,E] = S'^S'^ - S[S'^e 



I Ql -Aav'L 



0. 



(63) 
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The above expression can be solved numerically and the energy eigenvalues E could be 
obtained by setting Re[/[Q;, Vq, mo, mi, £■]] = and also Im[/[Q;, Vq, mo, mi, £■]] = since 
/[a, Vo, mo, mi, £■] = is complex. We search the numerical energy values for the interval 
m — Vo < E < m, since we are interested in bound states. Figs. 5 and 6 show the real energy 
eigenvalues for the PDM and constant mass cases, respectively. The eigenvalues which are 
shown with arrows are the points on the £^-axis where the Re[/[a, Vq, mo, mi, £']]-curve (solid 
line) and Im[/[Q;, T^o, mo, mi, £']]-curve (dotted line) cross. Prom the Figs. 5 and 6, one can 
see that the number of bound states in the case of PDM increases relative to the constant 
mass case. 

Finally, let us study our results in the case of low momentum limit which has become an 
attractive topic especially in the relativistic domain [3, 6, 14, 15]. In this limit, the Dirac 
equation has two distinct states: One has the energy E — m corresponding to particle state 
and the other one has the energy E — —m corresponding to anti-particle state where m is 
the particle mass [3]. 

Firstly, we investigate the resonance equation for the values of = ±m. For this case, 
Eq. (44) becomes 

-4at/L r(c)r(c -a-h) r(c2)r(c2 - 02 - &2) _ r(c)r(a + h-c) r(c2)r(a2 + &2 - C2) , . 

^ T{c-a)T{c-b)T{c2-a2)V{c2-h2)~ T{a)V{h) T{a2)T(h2) ' 

where the arguments are obtained from Eqs. (21) and (27) as 

a = v^E ±m) + cr ; 02 = i^{E ±m) + 5 , 

h = v{E ±m) - (7 + 1 ; 62 = t^{E ±m) -5 + 1, 

c = 1 ; C2 = 1 . (65) 

Secondly, we write the bound state equation given in Eq. (63) for E — ±m giving 
-WL r(cOr(c- -a'- h') r(4)r(4 -a',- h',) ^ T{c')T{a' + b'- d) nd,)V{a', + h',- 4) 

r(c' - a')T{c' - y) r(c'3 - a'^)T{d^ - fe'3) r(aOr(&o r(4)r(&'3) 

(66) 

where the arguments could be given from Eqs. (50) and (57) as 

a = u'{E ±m) + cr' ; 4 = u'{E ±m) + 6' , 

b' = u'{E ±m) -a' + l ; 63 = u'(E ±m) -5' + l, 

d^l- 4 = 1. (67) 
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Because of the substituting Vo — > —Vq, we see that v' — v^a' — a and 5' — 5, so 



a — a 



03 = 02 , 



b' ^ b ; 




c = c 



C3 = C2 . 



(68) 



which means that the resonance and bound state equations are equal for the low momentum 
limit in both of the constant and position-dependent mass cases. This result supports the 
one in Ref. [3] which declares that the conditions for the tunnelling without reflection 
trough a potential barrier V{x) of a Dirac particle with small momentum (resonances) and 
supporting of the potential well —V{x) a bound state energy E — ±m called supercriticality 
are the same. It should be noted that the transmission resonance of particles when they 
scatter off potential barriers is equivalent to the one of anti-particles when scattering off 
potential wells [6, 14, 15] and tunnelling of a Dirac particle trough a potential barrier is 
strongly related to the Klein paradox [14, 15]. 

V. CONCLUSIONS 

We have approximately solved the scattering and bound state problems in one- 
dimensional effective-mass Dirac equation for the Woods-Saxon potential and studied the 
problem by using an approximation in the mass distribution given as mi — >■ 0. Using this 
approximation, we have found reflection and transmission coefficients by analyzing the be- 
havior of the functions 0l(x) {x < 0) and (f)R{x) (x > 0) at x — )> =1=00 within the framework 
of position-dependent mass formalism. It has been observed that the coefficients oscillate 
for the range where E ^ mo and Vq » mo for both cases in which T and R change with 
energy and with potential parameter, respectively. The unitarity condition has also been 
checked numerically in the case of position-dependent mass. We have pointed out that the 
results for the case of constant mass arc similar with the ones obtained in the literature 
[3]. Meanwhile, the energy eigenvalue equation has been found by using the wave function 
obtained by imposing the boundary condition of a bound state. We have also studied the 
effect of the mass varying with coordinate on transmission resonances and on the results 
which are obtained for the low momentum limit (E — ±m). The relation between the trans- 
mission resonance and bound-state energy eigenvalues has been presented in the presence of 
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the position dependent mass case. 
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FIG. 1: transmission and reflection coeflicients versus energy with 
mo = 0.4, L = 10, a = 5, ^0 = 1-2. 
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FIG. 2: coefficient of transmission versus potential parameter Vq with 
mo = 0.4, L = 10,a = 5,E = 0.8. 




FIG. 3: transmission coefficient versus energy E with varying a (left plot: 
mo = 0.4, L = 10, Vo = 1.2) and varying L (right plot:mo = 0.4, a = 5, Vq = 1.2) in case of 

PDM (mi = 0.01). 
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FIG. 4: transmission coefficient versus potential parameter Vq with varying a (left plot: 
mo = 0.4, L = 10, E = 0.8) and varying L (right plot:mo = 0.4, a = 5, E = 0.8) in case of 

PDM (mi = 0.01). 
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FIG. 5: plots of Re[/[Q;, Vq, mo, mi, E]] and Im[/[Q;, V^), mo, mi, E]\ versus E for mi = 0.1 

(mo = 0.5, L = 5, q; = 10, Fo = 1). 
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